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Abstract  ^  We  introduce  the  fundamentals  of  cyclic 
multirate  systems  and  filter  banks  and  present  a  num¬ 
ber  of  important  differences  between  the  cyclic  and 
noncyclic  (traditional)  cases.  Some  of  the  additional 
freedom  offered  by  cyclic  systems  is  pointed  out,  and 
a  number  of  open  issues  are  summarized. 

1,  INTRODUCTION 

Digital  filter  banks  [l]-[2]  for  finite  length  signals  have 
been  considered  by  a  number  of  authors  [3]- [8].  Smith 
and  Eddins  [3]  introduced  this  idea  for  image  coding, 
and  among  other  significant  results,  proposed  replac¬ 
ing  linear  filtering  operations  with  cyclic  (or  circular) 
convolution.  They  also  noted  that  many  results  from 
the  noncyclic  case  carry  over  routinely.  In  this  paper 
we  draw  inspiration  from  these  references  to  develop 
the  theory  for  cyclic  filter  banks,  and  highlight  the  sig¬ 
nificant  differences  from  noncyclic  case.  A  cyclic (L) 
filter  bank  (with  the  letter  L  often  omitted)  has  an  in¬ 
put  x{n)  defined  only  for  0  <  n  <  L  —  1  (this  could  be 
a  symmetrically  extended  version  [3]  of  a  shorter  finite 
length  signal;  but  this  is  not  the  main  point  here). 

Fig.  1  shows  the  M-channel  cyclic(L)  filter  bank. 
We  assume  throughout  that 

L  =  KM,  K  =  integer.  (1) 

The  filters  have  impulse  responses  hi{n)  and  fi{n)  con¬ 
fined  toO<n<L  —  1,  and  their  I^point  DFTs 
are  denoted  as  Hi{k)  and  Fi{k).  All  convolutions  are 
cyclic(L).  So  we  can  regard  x{n)  to  be  periodic  (or 
cyclic)  with  period  L.  Since  a  cyclic(L)  signal  can  be 
represented  by  L  samples  of  the  Fourier  transform  (i.e., 
the  DFT  coefficients)  rather  than  the  entire  Fourier 
transform,  it  implies  more  freedom  in  the  design.  The 
paraunitary  and  power  complementary  properties,  and 
even  the  linear  phase  property  are  more  relaxed,  as 
they  are  imposed  only  on  a  discrete  frequency  grid. 

2.  BASICS  OF  CYCLIC  MULTIRATE  SYSTEMS 

The  notation  (or  just  W)  will  be  used 

throughout.  Since  the  frequency  variable  uj  is  re- 
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placed  by  the  discrete  version  27rkfL,  the  quantity 
Wj^  =  is  the  unit-delay  operator  analogous  to 

^-1  ^  The  cyclic  decimator,  denoted  by  J,  M  in 
Fig.  1,  has  the  input-output  relation  y{n)  =  x(Mn). 
With  x{n)  regarded  as  cyclic(L),  the  output  y{n)  is 
cyclic(if).  Thus  all  the  subband  signals  Vi{n)  are 
cyclic(K').  For  the  decimator,  let  X{k)  denote  the  L- 
point  DFT  of  x(n)  and  Y{k)  the  K'-point  DFT  of  its 
output  t/(n),  i.e., 

L-l 

A•(1b)  =  ^x(n)W^^  0<k<L-l  (2) 

n=0 

and  y(/k)  =  T,^~o  yin)W]^’",  0  <  <  AT  -  1.  It  can 

be  verified  that  Y(k)  =  X{Ki  +  k)/M,  for  0  < 

k  <  K  —  1.  The  cyclic  expander,  denoted  by  j  M  in 
Fig.  1,  has  a  periodic(i^)  input  x{n)  and  periodic(L) 
output  y{n)  related  by 

yin)  =  {  x(n/M),  n  =  mul.  of  M  (3) 
^  \  0  otherwise.  ^  ' 

The  corresponding  DFT  relation  is  Y{k)  =  X{k)  for 
0  <  A:  <  L  —  1.  Since  X{k)  are  the  AT-point  DFT  coef¬ 
ficients  of  x(n),  the  L-point  DFT  y(A:)  has  the  shorter 
period  K  =  L/M. 

Polyphase  representation.  The  L-point  DFT  of  a 
cyclic(L)  impulse  response  h{n)  can  be  expressed  as 

L-l  M-l 

H{k)  =  ^hin)Wf  =  J2^^L)'Ee{k)  (4) 

n=0  £=0 

where  Ee{k)  =  h{Mn  -h  for  0  <  A:  < 

L  —  1.  Thus,  Ei{k)  is  the  AT-point  DFT  of  the  ^th 

polyphase  component  e£(n)=h(Mn  .  From  the  def¬ 
inition  of  Ee{k)  we  see  that  it  is  cyclic(Ar).  Eq.  (4)  is 
analogous  to  the  traditional  Type  1  polyphase  decom¬ 
position  H{z)  =  z~^Ee{z^).  Fig.  2(b)  shows 

a  decimation  filter  redrawn  in  polyphase  form.  Simi¬ 
larly  the  Type  2  polyphase  form  is  given  by  H{k)  = 
W£^^Ri{k).  Since  the  polyphase  components 
Ei{k)  have  a  smaller  period  K,  they  can  be  relocated 
to  the  right  of  the  decimators  as  in  Fig.  2(c)  (similar 
to  the  use  of  noble  identities  in  the  noncyclic  case). 
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A  caution:  We  have  used  the  same  notation  Ei{k)  in 
Figs.  2(b)  and  2(c).  This  is  regarded  as  a  Jf-point 
DFT  in  Fig.  2(c)  and  an  2y-point  DFT  (with  values 
repeating  after  a  shorter  period  K)  in  Fig.  2(b). 

Return  now  to  the  analysis/synthesis  system  of 
Fig.  1.  With  the  filters  Hi{k)  represented  in  Type 
1  polyphase  form  and  the  filters  Fi{k)  in  Type  2  form, 
we  have  the  equivalent  representation  of  Fig.  3  where 
E(A:)  and  R(A:)  are  the  polyphase  matrices  of  the  cyclic 
filter  bank.  These  should  be  interpreted  as  RT-point 
DFTs,  e.g.,  E(/:)  =  12^0^  e{n)W^^,  The  filter  bank 
has  the  PR  property  iff  'R{k)E{k)  =  I  for  all  k. 

Orthonormality,  We  define  the  M-band  cyclic  fib 
ter  bank  to  be  orthonormal  if  the  matrix  E{k)  is 
unitary  for  all  k.  This  property  is  referred  to  as  the 
cyclic-paraunitary  property  The  perfect  reconstruc¬ 
tion  property  then  reduces  to  R(A;)  =  Et(A:),  or,  in 
terms  of  impulse  responses, 

fi{n)  =  h*{-n)  (5) 

(with  arguments  interpreted  modulo  L).  The  DFTs  are 
correspondingly  related  as  Fi{k)  =  H*{k),  Assuming 
x{n)  =  x{n)  in  Fig.  1  we  have 


M-lK-l 

^  ^  ^i(^)/i(^  0  <  n  <  L  —  1. 

^=0 

The  basis  functions  are  7]i^£{n)=fi{n  -  iM),  0  <  n  < 
L  —  1,  where  0<i<M  —  1,0<^<K  —  1.  Thus  there 
Sive  MK  ==  L  basis  functions  rii^e{n).  It  can  be  shown 
that  the  unitarity  of  R(/:)  is  equivalent  to  orthonor¬ 
mality  of  the  basis  This  orthonormality  can  be 

reexpressed  as  =  ^(^  ”  7n)6{£). 

As  in  traditional  filter  banks,  orthonormality  of  the 
cyclic  filter  bank  implies  the  unit  energy  property 
J2n=o  l/i(^)P  =  1  and  the  power  complementary  prop- 
ertyEiloMi^iWP  =  M. 

3.  CYCLIC  VERSUS  NONCYCLIC  FILTER  BANKS 

In  the  cyclic(L)  case,  any  transfer  function  can  be  ex¬ 
pressed  in  the  form  H{k)  =  X]n=o  ^(^)^l”-  The  non- 
cyclic  counterpart  of  H{k)  is  defined  as  Hnon{z)  = 
^  interpolated  version  in  the 
frequency  domain,  since  H{k)  =  samples  of  Hnon{z) 
aX  z  =  Similarly  the  noncyclic  counterpart 

of  E{k)  =  El"o  e(n)W^"  (Fig.  3)  is  Enon{z)  = 
J2n=o  e(n)2:“".  If  Enon(z)  IS  paraunitary  (PU)  it  read- 
ily  follows  that  E(fc)  is  cyclic-PU  because  each  k  cor¬ 
responds  to  a  special  2:  on  the  unit  circle.  But  the 
converse  does  not  hold  as  we  shall  see.  Thus  the 
cyclic-PU  property  is  less  stringent  than  traditional 
PU.  More  generally,  if  we  impose  a  certain  constraint 
on  the  cyclic  system,  it  does  not  mean  that  the  non¬ 
cyclic  counterpart  has  to  satisfy  the  same  constraint. 


To  demonstrate,  consider  the  second  order  cyclic(4) 
transfer  function 


G(k)  -  0.5  +  0.5(i  -  1) Wf  +  0.5i  Wf 

Using  the  facts  that  =  1  and  =  —1,  it  is  readily 
verified  that  |G(A:)p  =  1  for  all  k.  Thus  G{k)  is  allpass 
in  the  cyclic(4)  sense.  However,  the  noncyclic  version 
Gnon{z)  =  0.5  -h  0.5{j  *“  1)2;“^  -h  0.bjz~^  which  is  an 
FIR  filter,  is  evidently  not  allpass.  If  we  now  construct 
the  M  X  M  polyphase  matrix  E{k)  =  G{k)lM  (for 
arbitrary  M)  then  it  is  cyclic-PU,  but  the  noncyclic 
version  Enon(2)  =  Gnon{z)I  is  not  PU.  As  a  second 
example,  consider  the  cyclic(3)  analysis  bank 


H(fc)  = 


fi/oWl 

hM 

mk). 


=  ao  + 311^3''+ 32 1^2'= 


where  are  the  column- vectors  given  by 


1 

T 

1 

■  1  ■ 

1 

T 

Vio 

1 

0 

--1 

-2 

11 

(N 

0 

1 

It  can  be  verified  that  H^{k)ll{k)  =  1,  that  is,  the 
three  transfer  functions  Ho{k),  Hi{k)  and  i?2(^)  satisfy 
the  cyclic  power  complementary  property  \Ho{k)\'^  -{- 
|^^i(^)P  +  \H2{k)\'^  =  1.  But  the  noncyclic  version 
Hnon(2)  =  power  Complementary, 

that  is,  Hnon(2)H 

non{^)  ^  !• 

Cyclic  allpass  filters,  A  cyclic(L)  allpass  filter  satis¬ 
fies  H{k)  =  for  0  <  A:  <  L  —  1.  We  can  always 

write  Hik)  =  c(Elo«’N-nW^^)/(Elo«'nW2")  with 
\c\  =  1.  For  example,  we  can  let  A”  =  L  —  1  and  set 
e-i0(fc)/2  _  The  coefficients  bn  are  es¬ 
sentially  the  inverse  DFT  coefficients  of  and 

can  readily  be  identified.  More  interesting  is  the  open 
problem  of  obtaining  the  rational  form  with  smallest 
order  N.  Now  consider  the  L  x  L  circulant  matrix  H 
whose  first  row  is  /i(0),  /i(l),  h{2) . . .  We  know  that  cir¬ 
culant  matrices  are  diagonalized  by  the  DFT  matrix 
(which  is  unitary)  and  that  the  eigenvalues  are  the 
DFT  coefficients  H{k),  Using  this  we  see  that  the  all¬ 
pass  property  is  equivalent  to  unitariness  o/H. 

Factorization  of  cyclic  paraunitary  systems,  Non¬ 
cyclic  PU  systems  can  be  factored  [2]  in  terms  of  the 

building  blocks  I  —  u^uj  -f  2:“^Uiut,  where  utu^  = 
1.  But  in  the  cyclic  case,  factorization  in  terms  of 

(I-u^ut  -hVF^Uiut)  is  not  always  possible.  When  it  is 
possible,  we  can  replace  with  z~^  in  the  factoriza¬ 
tion,  and  obtain  a  noncyclic  PU  counterpart  Enon{z). 

Determinant  of  cyclic  paraunitary  system.  If  E(A:)  is 
cyclic  PU,  then  [det  E(A:)]  =  which  is  allpass. 

In  the  noncyclic  FIR  case  the  degree  of  [det  Enon{z)] 


2 


is  equal  to  the  degree  of  Enon{z),  but  the  same  is  not 
true  in  the  cyclic  case.  Thus  let 


E{k)^ 


cos  0{k) 
-sin  0{k) 


sm6(^k) 
cosOlk)  ’ 


0  <  A:  <  L  -  1  (6) 


which  is  cyclic  PU  with  [det  E{k)]  =  1.  The  determi- 
nantal  degree  is  zero,  regardless  of  the  degree  of  E(fc) 
(which  is  the  minimum  number  of  unit  delay  elements 
required  in  the  implementation). 

Nyquist  and  Linear-Phase  Properties.  Consider  a 
cyclic(6)  transfer  function 


that  g{n)  =  0  for  iV  <  n  <  L  -  A^)  the  cyclic  spec¬ 
tral  factor  HQ{k)  may  not  be  FIR  (i.e.,  ho{k)  could  be 
nonzero  for  all  n).  See  Fig.  4.  An  open  question  here 
is,  what  is  the  most  general  form  of  the  phase  response 
<i>{k)  which  ensures  that  Ho{k)  is  FIR  (possibly  with 
length  AT  -fl). 

4.  CYCLIC  DIFFERENCE  EQUATIONS 

Consider  a  cyclic(L)  transfer  function  of  the  form 
H{k)  =  (oo  -f  aiW^)/(l  -  hWl),  The  input  and  out¬ 
put  of  this  system  are  constrained  by  {l  —  bW^)Y{k)  = 
(ao  +  aiWj^)X{k),  By  taking  inverse  DFT  we  obtain 


Hoik)  =  (l  +  -  Wi'^  +  Wt  +  Wt)/Vs. 

We  can  regard  this  as  FIR  in  the  sense  that  h{n) 
is  nonzero  only  on  a  subset  of  points  in  0  <  n  < 
L  —  1.  (In  the  cyclic  case  this  is  the  only  FIR-definition 
that  makes  sense).  The  symmetry  of  the  impulse 
response  implies  the  linear-phase  property,  that  is, 
Ho{k)  =  X  R{k)  where  R{k)  is  real.  Now  consider 
G{k)  =  |Fo(A:)p.  Its  cyclic(6)  impulse  response  g{n) 
satisfies  the  halfband  property  g{2n)  =  6{n),  Equiva¬ 
lently  H{k)  is  power  symmetric  [2]  in  the  cyclic  sense: 

\Hoik)f  +  \Ho{i  +  k)\^  ^2,  0<fc<5.  (7) 

So  we  have  a  cyclic(6)  filter  Ho{k)  which  is  both  linear- 
phase  and  power-symmetric.  This  is  not  possible  in  the 
noncyclic  FIR  case  [2]. 

Using  the  above  example  we  can  construct  a  two- 
channel  cyclic (6)  orthonormal  filter  bank  where  the  fil¬ 
ters  are  nontrivial  linear  phase  filters.  (Such  construc¬ 
tions  are  not  possible  in  the  noncyclic  FIR  case  [2].) 
For  this  choose  ho{n)  as  above,  and  the  remaining  three 
filters  as  hi{n)  =  (— l)"^/io(l  -  n),  fi{n)  =  /j*(-n). 
That  is,  Hi{k)  =  -W^H^{k-K)  where  K  =  L/2,  and 
Fi{k)  =  H*{k)  (recall  all  arguments  are  interpreted 
modulo-6).  This  an  example  of  a  cyclic(6)  version  of 
the  CQF  design. 

More  on  the  CQF  Design.  Given  a  cyclic(L)  half¬ 
band  filter  g[n)  with  the  property  G{k)  >  0,  we  can 
construct  infinitely  many  cyclic(L)  two  channel  or¬ 
thonormal  filter  banks  by  choosing  Ho{k)  as  a  spectral 
factor 

=  (8) 

and  the  remaining  three  filters  according  to  the  CQF 
equations  above.  Recall  that  for  noncyclic  FIR  fil¬ 
ters,  the  usual  definition  of  spectral  factors  allows  only 
finitely  many  phase  responses,  and  the  spectral  factors 
all  have  the  same  length  AT -1-1.  In  Eq.  (8)  however,  the 
phase  response  4){k)  of  the  spectral  factor  Ho{k)  is  ar¬ 
bitrary.  In  particular,  the  choice  4>{k)  =  0  would  yield 
linear  phase  analysis  filters  Ho{k)  and  Hi{k)  with  good 
frequency  responses  (if  G{k)  is  a  good  lowpass  filter). 
However,  even  if  G{k)  is  cyclic(L)  FIR  (in  the  sense 


y(n)  =  by{n  -  1)  +  aox{n)  -f  aix{n  -  1)  (9) 

Since  the  time-indices  are  interpreted  modulo  L,  this 
is  a  cyclic  difference  equation.  To  demonstrate  what 
is  involved,  let  L  =  3.  Repeated  use  of  the  d.e.  (and 
the  facts  that  y{2)  =  2/(-l),  2/(3)  =  y(0)  etc.)  there¬ 
fore  yields  the  conclusion  that  (1  -  &^)2/(0)  =  aoa;(0)  -{- 
aix{l)-\-a2x{2)  for  some  constants  ai.  Thus  the  initial 
condition  y(0)  is  not  arbitrary,  but  is  uniquely  deter¬ 
mined  as  long  as  6^  7^  1,  that  is,  as  long  as  6  7^ 
for  any  i.  This  is  equivalent  to  the  obvious  require¬ 
ment  that  the  denominator  in  H{k)  does  not  become 
zero  for  any  k.  We  can  then  write  Y{k)  =  H{k)X{k), 
and  the  inverse  DFT  y{n)  is  determined  for  all  n. 
More  generally,  consider  the  cyclic(L)  transfer  func- 

tion  Hik)  =  (Elo«nW^r)/(l  +  Eli  bnWtn-  With 
the  implicit  assumption  that  the  denominator  does  not 
vanish  for  any  k,  the  output  is  fully  determined  by  the 
input.  In  particular  the  ^%nitial”  condition  is  predeter¬ 
mined,  and  is  not  arbitrary. 

To  compute  y(n)  efficiently,  we  can  use  a  recursive 
d.e.  similar  to  (9)  provided  we  know  how  to  identify 
the  (unique)  initial  condition.  A  convenient  way  is  the 
state  space  apprach.  From  the  direct-form  structure  of 
Fig.  5  we  can  identify  a  set  of  N  state  variables  Vi{n) 
(outputs  of  the  unit  delay  elements  W^)  and  obtain 
equations  of  the  form 

v(n  -h  1)  =  Av(n)  +  Bx(n)  (10) 

and  y{n)  =  Cv(n)  -f  Dx{n)  where  v(n)  is  the  state 
vector  with  components  Vi{n).  Repeated  use  of  (10) 
yields  v(L)  =  A^v(O)  -1-  a  linear  combination  of  sam¬ 
ples  of  x{n).  Since  all  the  time-indices  are  interpreted 
modulo-L,  we  have  v(L)  =  v(0),  so  (I  -  A^)v(O)  = 
linear  combination  of  samples  of  x(n).  Thus  we  can 
identify  the  initial  state  v(0),  if  I  -  A^  is  nonsingular 
(i.e.,  no  eigenvalue  of  A  has  the  form  Wl  for  any  i). 

As  in  the  noncyclic  case  the  transfer  function  H{k) 
can  be  expressed  as  H{k)  =  D-\-C{W£^l—A)"^By  but 
the  impulse  response  is  h{n)  =  CA’^~^(I  — A^)“^B  for 

0  <  n  <  L  and  h(0)  =  D+ (^“ ^ 

is  because  the  initial  condition  v(0)  is  predetermined 
and  cannot  be  set  to  zero. 
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The  system  matrix  for  an  LTI  cyclic  implementa¬ 


tion  is  defined  as 


A  B 
C  D 


.  If  this  is  unitary,  one 


can  verify  (as  in  the  noncyclic  case  [2])  that  E{k)  is 
cyclic  PU.  However  the  converse  is  not  true  in  the 
cyclic  case:  even  if  E(A:)  is  PU,  there  may  not  ex¬ 
ist  a  minimal  structure  with  unitary  system  matrix. 
Whenever  it  does,  we  can  find  an  interpolated  PU  ma¬ 
trix  Enon{z)  by  replacing  Wj^  by  in  the  expression 
H{k)  =  D  +  C{W£^I  -  A)“^B,  and  then  factorize  it. 
In  the  factored  form  if  we  replace  z~^  with  we 
obtain  a  factorization  of  E(A;). 


Fig.  1 .  The  cyclic  filter  bank. 


5.  CONCLUDING  REMARKS 

Perhaps  the  question  of  greatest  interest  for  future 
work  is,  how  to  exploit  the  extra  freedom  offered  by 
the  cyclic  system,  in  the  design  of  subband  coders. 
Can  this  be  exploited  to  obtain  increased  coding  gain 
(or  compression),  or  to  obtain  reduced  complexity  of 
implementation?  These  require  further  study.  The  no 
tion  of  cyclic  LTI  systems  also  opens  up  other  problems 
in  the  general  arena  of  signal  and  system  theory. 
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Fig.  4.  A  cycIic(L)  FIR  filter  with  non-FlR  spectral  factor. 


Fig.  2.  (a)  A  decimation  filter  and  {b),(c)  polyphase  forms. 


Fig.  3.  Polyphase  form  of  the  cyclic  filter  bank. 


x(n) 


The  system  matrix  for  an  LTI  cyclic  implementa¬ 


tion  is  defined  as 


A  B 
C  D 


.  If  this  is  unitary,  one 


can  verify  (as  in  the  noncyclic  case  [2])  that  E(fc)  is 
cyclic  PU.  However  the  converse  is  not  true  in  the 
cyclic  case:  even  if  E{k)  is  PU,  there  may  not  ex¬ 
ist  a  minimal  structure  with  unitary  system  matrix. 
Whenever  it  does,  we  can  find  an  interpolated  PU  ma¬ 
trix  Enon{z)  by  replacing  Wj^  by  in  the  expression 
H{k)  =  D  +  C{W^^1  -  and  then  factorize  it. 

In  the  factored  form  if  we  replace  with  we 
obtain  a  factorization  of  E{k). 


Rg.  1 .  The  cydic  filter  bank. 


5.  CONCLUDING  REMARKS 

Perhaps  the  question  of  greatest  interest  for  future 
work  is,  how  to  exploit  the  extra  freedom  offered  by 
the  cyclic  system,  in  the  design  of  subband  coders. 
Can  this  be  exploited  to  obtain  increased  coding  gain 
(or  compression),  or  to  obtain  reduced  complexity  of 
implementation?  These  require  further  study.  The  no¬ 
tion  of  cyclic  LTI  systems  also  opens  up  other  problems 
in  the  general  arena  of  signal  and  system  theory. 
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